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We divide the circular boundary of a hyperbolic lattice into four equal intervals, and study the 
probability of a percolation crossing between an opposite pair, as a function of the bond occupation 
probability p. We consider the {7,3} (heptagonal), enhanced or extended binary tree (EBT), the 
EBT-dual, and {5,5} (pentagonal) lattices. We find that the crossing probability increases gradually 
from zero to one as p increases from the lower pi to the upper Pu critical values. We find bounds 
and estimates for the values of pi and pu for these lattices, and identify the self-duality point 
p* corresponding to where the crossing probability equals 1/2. Comparison is made with recent 
numerical and theoretical results. 

PACS numbers: 64.60.ah, 64.60.De, 05.50.+q 



I. INTRODUCTION 

Hyperbolic lattices represent curved surfaces in a space 
that is effectively of infinite dimensions. While long of 
interest to mathematicians [1 , and even artists [2 , such 
lattices have only relatively recently been studied in sta- 
tistical physics, where many problems |3Hl3j including 
percolation ^141427] have been examined. Investigation 
has also been carried out on closely-related hierarchical 
lattices [28ti3T] . The study of hyperbolic lattices helps in 
the understanding of how geometry affects the behavior 
of systems. There are also physical systems that show 
negative curvature on the nanoscale [32 . Networks are 
not confined to a physical dimensionality and can show 
hyperbolic behavior, and the current strong interest in 
network physics is another motivation to study these sys- 
tems. 

Connectivity in networks is described by percolation, 
which has been studied for a wide variety of systems for 
well over 50 years [33]. For most systems in percolation, 
there is typically a bond or site occupation probability 

such that when p is less than a threshold value all 
connected components are finite, while above Pc there is 
an infinite network which connects every region of the 
system. 

For hyperbolic lattices, however, percolation shows a 
more complicated behavior, with two distinct transitions, 
as described below. Various numerical and theoretical 
methods have been developed to study these transitions, 
and in particular to find the threshold values [14-27J. 

In this paper, we investigate percolation on hyperbolic 
lattices by considering the crossing probability, a tech- 
nique which has not been applied to this system before. 
The lattices we study are some that have been consid- 
ered by others, so that a comparison of the values for 
the transition points can be made. We also study (for 
the first time numerically) percolation on the pentagonal 
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lattice, which has some interesting self-dual features. In 
general, the determination of the transition points allows 
one to compare the different numerical methods and to 
test theoretical predictions and bounds. These threshold 
points will be useful for future studies of the nature of the 
transition behavior, such as the determination of critical 
exponents. 




FIG. 1. (color online) Heptangonal lattice {7,3} (black or 
dark) and the dual lattice {3, 7} (magenta or light). 

Common examples of hyperbolic lattices are those 
composed of identical polygons of n sides, m of which 
meet at a vertex. These lattices can be character- 
ized by the Schlafli symbol {n,m}, corresponding to the 
Griinbaum-Shepard [34] notation (n^). Thus, {6,3} is 
a regular planar hexagonal (honeycomb) lattice, while 
{7, 3} is the heptagonal hyperbolic lattice shown in Fig. 
[l] The dual to the heptagonal lattice is the {3, 7} lattice, 
also shown in Fig. [l] The self-dual hyperbolic {5,5} is 
shown in Fig. [9| 

Recently, another type of hyperbolic lattice has been 
introduced, the Enhanced or Extended Binary Tree 
(EBT) [14| [T8|), which is made by adding transverse 
bonds to the Bethe lattice. The EBT, which is simpler 
to represent and code on a computer, has been studied 
extensively for percolation [T8H2Q[ l22ti24] . The EBT and 
its dual are shown in Fig. |2j 
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FIG. 2. (color online) The Enhanced Binary Tree (EBT) lat- 
tice (black or dark) and EBT-dual lattice (red or light). The 
EBT is made by joining four trees together. We joined four 
bonds at the center rather than having two join there as in 
[241. 



The general picture that has emerged for percolation 
on hyperbolic lattices p31 HSl [13 HH] is that there are 
two thresholds pi and Pu , and the behavior is continuous 
between them. For < p < Pu there are no "infinite" 
(large) clusters connecting the central area to the bound- 
ary sites, for the intermediate region pi < p < Pu there 
are many infinite clusters touching the boundary, and for 
Pu < p < 1, there is exactly one infinite cluster. These 
three regions persist in the limit that the system size goes 
to infinity. This behavior is in contrast to ordinary per- 
colation, in which there is no intermediate region (in an 
infinite system) and the crossing behavior is discontinu- 
ous. 

Here we study percolation on hyperbolic lattices by ex- 
amining a suitably defined crossing probability R{p) as a 
function of the bond occupation probability p. The cross- 
ing probability is often studied in ordinary percolation to 
locate the threshold and to investigate the critical scal- 
ing behavior [35 -37 . There, the crossing probability (the 
probability of a continuous path from one opposite side to 
another) becomes steeper as the size of the system is in- 
creased, with a slope proportional to L^/^, where u is the 
correlation- length exponent, equal to 4/3 in two dimen- 
sions. This behavior defines the transition point uniquely 
in the limit of L ^ oo [33]. Furthermore, when the sys- 
tem boundary is symmetric such as a perfect square, the 
crossing probability between opposite sides is exactly 1/2 
(because the dual lattice percolates if the original lattice 
does not) [38l|39]. A crossing probability of 1/2 applies 
to a disk also, with the circumference divided into four 
equal intervals. In this paper, we set up a similar crossing 
problem for hyperbolic system by dividing the bound- 
ary of a finite system into four equal-size intervals (or as 
equal as possible), and study the probability of crossing 
between an opposite pair of these intervals. We consider 
the heptagonal, EBT, EBT-dual, and pentagonal hyper- 
bolic lattices, and investigate how the resulting crossing 
probability behaves with p. We also discuss how the value 



of p = p* that corresponds to a crossing probability of 
exactly 1/2 relates to the two transition points pi and Pu- 
In the following sections we discuss the methods (II), 
the results (III), and the conclusions (IV). 



II. METHOD 

We begin by generating a hyperbolic lattice to a fixed 
number of generations or levels, so that the outside is 
essentially circular as in Fig. [l] Practically, it is only 
feasible to generate a relatively small number of levels 
(up to 10 - 15) because of the exponential growth in the 
number of lattice sites with level. 

For the heptagonal lattice with an open heptagon cen- 
tered at the origin, we can derive an expression for the 
number of sites N{1) as a function of level / as follows: 
Let ai equal the number of new sites which connect to 
the next generation, and bi equal the number of new sites 
which connect to the previous generation. Inspection of 
the Fig. [l] shows that we have the relations 



ai 
hi 



3a^_i 



(1) 



Thus, the total number of sites up to level / is equal to 
Yl\'=i{^i + ^0- means of generating function tech- 
niques, we find the explicit relation 



N{1) = 7 



^3 + ^5' 



(2) 



which yields N{1) = 7, 35, 112, 315, 847, 2240, 5887, 
15 435, 40 432, 105 875. . . for / = 1, 2, . . . , 10, . . .. In Ref. 
[17], the corresponding formula for N{1) with a vertex 
rather than an open heptagon at the center of the system 
is given. For large N{1) grows exponentially as ~ 7[(3-h 
^/b)/2y. These N{1) are related to other mathematical 
quantities, such as the number of fixed points of period / 
in iterations of Arnold's cat map at its hyperbolic fixed 
point, multiplied by 7 [40 . 

For the pentagonal lattice, we find 



ai = 5a/_i 
bi 3a/_i 



^bi 
2bi 



which yields 



N{1) 



(3) 



(4) 



and equals 5, 45, 320, 2205, 15125, 103 680, 710 645, 
4 870 845. . . for / = 1, 2, . . . , 8, . . .. Here the N{1) are re- 
lated to the Fibonacci numbers F{1) by N{1) = 5-F(2/)^ 
For the EBT, we consider a geometry with four trees 
meeting at the origin as shown in Fig. [2| so that it is 
easy to divide the boundary into four equal intervals. 



3 



The number of sites grows as A^(/) = 2^+^ — 3. For the 
EBT-dual, we have N{1) = -2^+2 - 4. 

We apphed the algorithm of [411 HZ^ to find the cross- 
ing probabihty on these four lattices. This algorithm 
allows one to find an estimate of R{p) for all values of 
p in a single sweep of the lattice; averaging over many 
sweeps yields an accurate estimate of R{p). The connec- 
tions between points in a cluster are represented by a tree 
structure, bonds are added one at a time, and clusters are 
joined together by means of a union-find operation. The 
algorithm yields the crossing probability Rn as a function 
of the number of occupied bonds n added to the system, 
corresponding to a fixed-n or canonical ensemble. To get 
the grand canonical result R{p) corresponding to a fixed 
probability p, one must convolve Rn with the binomial 
distribution: 

where N is the total number of bonds in the system. 
For large and the binomial distribution becomes very 
sharp, and for many problems it is not necessary to carry 
out this convolution, but instead use just the value at the 
maximum of the distribution n = Np, so that R{p) ^ 
Rnp' However, for smaller systems it is necessary to use 
this convolution to get accurate results. 

III. RESULTS 

We carried out simulations for each of the four lattices, 
recording R{p) at 500 values of p. Below we describe the 
results for each lattice. 



A. Heptagonal {7,3} lattice 

We considered the heptagonal {7, 3} lattice up to level 
/ = 10 with N{1) = 105 875 total sites. Fig. [3] shows 
the resulting R{p) as a function of p for levels 5, . . . , 10. 
We find a gradual transition of R{p) from to 1 as the 
p increases, as is typical for finite systems for ordinary 
percolation. However, here the width of the transition 
region is more spread out and, more significantly, the 
width limits to a non-zero value as / -> oo. Equivalently 
the slope at the infiection point limits to a finite value 
as / ^ oo. In Fig. |4] we plot the maximum slope as 
a function of A/'(/)~^ ' where N{1) is given by ([2|, and 
find an extrapolation to a maximum value of ~ 6.12. 
The exponent —0.7 was chosen empirically to get a fairly 
straight line; different choices of the exponent do not 
change the intercept significantly and especially do not 
change the conclusion that the slope limits to a finite 
value as ^ oo. 

The close-up in Fig. [3] shows that the curves do not 
quite cross at a single point, but the crossing point 
changes with N. For a perfectly self-dual system in which 
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FIG. 3. (color online) Curves of crossing probability for the 
heptagonal lattice, convoluted with ([s]), for various levels /; 
the slope increases as / increases. The dashed line passes 
through the inflection point, and its intercepts with the lines 
at R — and R — 1 gives our estimates of pf and - Inset: 
close-up near crossing point. 
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FIG. 4. Maximum of the slope of R{p) versus N for the 
heptagonal lattice, where N is given by (pi). 



the dual lattice is identical to the original lattice, such as 
bond percolation on a square lattice and square bound- 
ary in ordinary percolation, the curves cross at a single 
point corresponding to R = 1/2 and p = 1/2, but be- 
cause this system is not self-dual, one would not expect 
the crossing to be at (1/2, 1/2) here. 

We define the duality point (/) as the value of p where 
R{p) = 1/2. We call this the duality point because on 
a truly dual lattice the crossing probability should also 
be 1/2 (although below we see that there are differences 
in the center that limit the extent that one can make a 
completely self-dual system). We find p*(oo) ~ 0.6759 
by extrapolating to / = oo as shown in Fig. [5] Here we 
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FIG. 5. (color online) Dual point p*{l) versus 1/A^(/), where 
N(l) is the number of sites on the lattices of levels / = 5, . . . , 10 
for the heptagonal lattice as given by We show results 
for the raw (canonical) (upper points) and convoluted (grand 
canonical) (lower points) data; both extrapolate to the same 
value, p ^ 0.6759, as L ^ oo. 

observed a scaling of order 1/A^(/). 

The transition points for the heptagonal lattice were 
found by Back et al. [17 to be pi ^ 0.53 and Pu ^ 0.72, 
and on the dual lattice {3,7} they found pi ~ 0.20 and 
Pu ^ 0.37. These four values are evidently not completely 
consistent because one should have, for any lattice and 
its dual, 

Pn+Pr^ = l (6) 

One would expect that for p > Pu^ R{p) = 1, and 
for p < Pu R{p) = 0- However, how R{p) approaches 
those values from the region pi < p < pu is not clear. It 
appears from our data the approach is tangential (with 
slope zero), and therefore it is rather hard to identify 
the transition points accurately. We can, however, find 
bounds to that behavior by drawing a tangent line from 
the inflection point. Drawing a line through the inflection 
point in Fig. [3] with the maximum slope ~ 6.12, and the 
intercepts for R{p) = and R{p) = 1 give us the rather 
crude bounds Pi < pf = 0.594 and Pu > Vu — 0-758. 

A more precise method to get bounds or estimates for 
the transition point is to look at the values of p where 
R{p) = e and R{p) = 1 — e, where we chose e = 10~^, and 
then extrapolating to L = oo. Fig. [6] shows that these 
estimates appear to scale as 1//, and extrapolating the 
points to / ^ oo gives the values of and listed in 
Table [T[ In principle, the smaller value of e the better, 
but noise in the data and precision of the numbers in our 
output flies limited how small we could make e. The flxed 
number of points that we recorded (500) also limited the 
flnal precision of the thresholds. Varying the value of e 
by an order of magnitude in each direction and flnding 



the extrapolated thresholds, we estimate that the overall 
error in our threshold estimates is about ±0.01. 
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FIG. 6. Values of p where R{p) — 1 — 10~^ (upper data points) 
and R(p) — 10~^ (lower data points), plotted as a function 
1//, for the {7,3} heptagonal lattice. The linear extrapolation 
to / ^ oo gives our estimates p^ — 0.810 and pi — 0.551. Ex- 
trapolations for the other lattices show similar linear behavior 
and the values for and p J are given in Table [l] 



B. EBT and EBT-dual lattices 

We simulated the EBT lattice to the level of 15, and 
the EBT-dual lattice to the level of 10. Figs. [7| and [8] 
show the resulting crossing probability distribution for 
these two lattices. For the EBT, the maximum slope 
converges to ^ 6.79. Its duality point is at p* ^ 0.4299, 
yielding the bounds pf ~ 0.356 and p^ ~ 0.503. The 
EBT-dual's crossing probability distribution curve also 
converges to a maximum slope 6.83, the duality point 
of which is at p* ^ 0.5698, yielding the bounds pf = 
0.497 and p^ = 0.643. These bounds satisfy the expected 
duality ([6| within errors. The estimates are also found to 
the scale as 1// and the resulting values pf and given 
in Table [11 These estimates do not satisfy the duality 
relations l6| very precisely, reflecting rather large error 
bars in their values. 



C. Pentagonal lattice 

We also considered the pentagonal {5, 5} lattice, which 
is shown in Fig. [9j This lattice is interesting because it 
is self-dual in an inflnite system. For the systems of a 
flnite number of levels, it is not precisely self-dual be- 
cause the center is different: on what we call the pentag- 
onal lattice, there is a pentagon at the center, while for 
the pentagonal-dual, there is a vertex at the center (see 
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FIG. 9. (color online) Pentagonal (black or dark) and dual 
pentagonal (red or light), both {5,5}. 



FIG. 7. (color online) The crossing probability as a function 
of p for the EBT lattice, for systems of 5 - 15 levels. The 
dashed line passes through the inflection point. 
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FIG. 10. (color online) The crossing probability as a func- 
tion of p for the {5,5} or pentagonal lattice, for systems of 
5, 6, and 7 levels. The curves are nearly indistinguishable on 
this plot. The dashed line passes through the inflection point. 



FIG. 8. (color online) The crossing probability as a function 
of p for the EBT-dual lattice, for systems of 5 - 10 levels. The 
dashed line passes through the inflection point. 



Fig[lO|. We find that p* = 0.506 ± 0.001, so it is not ex- 
actly at 0.5 as one might expect from duality. Evidently, 
the central region plays an important role and a signif- 
icant fraction of the percolating clusters connecting op- 
posite sides pass through it, making the pentagonal and 
pentagonal-dual lattices slightly different with respect to 
the crossing problem we consider. 

The slope of its crossing probability curve converges to 
a maximum value ~ 3.12, with bounds pf ~ 0.346 and 
p^ ^ 0.666, which indicates the distribution is nearly 
symmetric. 

Recently, Delfosse and Zemor [43] have shown that, for 
any self-dual hyperbolic lattice {m, m}, l/(m — 1) < pi < 
2/m, so that for m = 5, 1/A<pi < 2/5. Our bounds pf 



and p^ fall well within these values, and our estimates p^ 
and are close to the bound 1/4 and (by duality) 3/4, 
respectively. This bound follows from approximating the 
hyperbolic lattice as a tree (Bethe lattice) of coordination 
number 5. 



Lattice 


* 

P 


max. slope 




Pu 


Pt 


pI 


{7,3} 


0.6759 


6.12 


0.594 


0.758 


0.551 


0.810 


EBT 


0.4299 


6.79 


0.356 


0.503 


0.306 


0.564 


EBT-dual 


0.5698 


6.83 


0.497 


0.643 


0.452 


0.699 


{5,5} 


0.506 


3.12 


0.346 


0.666 


0.263 


0.749 



TABLE I. Values of the dual point p* where R{p) = 1/2, the 
slope at that point, our bounds pf and p^ for the various 
lattices we studied, such that pi < pf , and pu > pf , and our 
extrapolated estimates of the transition points pf and p^ ■ In 
general, the numbers are expected to be accurate to about ±1 
in the last digit shown, except the estimates and pj, which 
are expected to be accurate to about ±10 in the last digits. 
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Lattice 


Pi 


Pu 


Ref. 


{7,3} 


0.53 


0.72 


LlZi 


{3,7} 


0.20 


0.37 


m 


EBT 


0.304(1) 


0.564(1) 


CH] 


77 




0.48 


m 




(7l3 -3)/2 ^ 0.3028 


0.5 


[23^ 


EBT-dual 


0.436(1) 


0.696(1) 




{5,5} 


0.25 <pi< 0.4 




|43| 



TABLE IL Previous values of the transition points. 



IV. CONCLUSIONS 

In summary, we have the following results and conclu- 
sions: 

(i) The crossing probability approaches a continuous 
S-shaped curve with a finite maximum slope at the in- 
flection point as / ^ oo. 

(ii) By drawing a tangent line through the inflection 
point and finding its intercept with R{p) = and R{p) = 
1, we find the bounds pf and p^ for the transition points 
pi and Pu listed in Table [l[ Also, by extrapolating where 
R{p) — e and R{p) = 1 — e to L ^ oo, we find the 
estimates pf and p^, also listed in Table [l| In comparison, 
previously measured and predicted values of pi and Pu 
(determined through other methods) are listed in Table 

nn 

(iii) For the {7,3} lattice, the reported value Pu = 0.72 
[17] is inconsistent with our lower bound p^ = 0.758 
and estimate 0.810. However, those authors' value for 

= 0.20 on the dual lattice {3,7} is consistent with this 
bound, by ([6|. 

(iv) For the EBT lattice, our bound p^ = 0.503 and 
especially our estimate p^ = 0.564 are inconsistent with 
the prediction pu = 1/2 [23 . Our results for pu are 
however consistent with the measurement of Pu by [18 . 
For pi^ our estimate 0.306 is in substantial agreement 
with the results of both Refs. [17 and [18 . 

(v) For the EBT-dual lattice, our bounds and estimates 
for the transition points agree with those of [18] within 
expected errors. 

(vi) For the {5,5} lattice, we report measurements of 
the thresholds for the first time, and our estimates are 
close to the theoretical bounds pi = 1/4 and Pu = 3/4, 



which follow by assuming a tree structure [43] . 

(vii) We determine the point where i?(p*) = 1/2 for 
all four lattices we consider, and find that the behavior of 
R{p) is nearly symmetric about that point. For the {5,5} 
lattice, p* ~ 0.506, slightly larger than the value 0.5 one 
might expect from self-duality. We believe the deviation 
from 0.5 is due to the non-equivalent configurations at 
the center for the lattice and its dual, implying that the 
two finite systems are not exactly dual to each other. 

(viii) We have found a variety of finite-size scaling re- 
lations, as given in Figs. [4] |5| and [6| Another scaling 
relation that exists in the literature is that of Nogawa 
and Hasegawa [18 , who find that the mass of the root 
cluster in the EBT scales as N'^^p^ in the intermediate 
region, where is a function of p. Clearly it would be 
desirable to have a general scaling theory that combines 
all of these finite-size scaling relations. This is an inter- 
esting problem for future study. 

(ix) Another area for future study is site percola- 
tion on hyperbolic lattices. For site percolation on fully 
triangulated lattices in ordinary two-dimensional space, 
Pc = 1/2. For fully triangulated hyperbolic lattices, such 
as the {3, 7} lattice, we guess that the behavior of R{p) 
will be precisely symmetric about p = 1/2, because of its 
self-matching property. Also, because bond percolation 
on a given lattice is equivalent to site percolation on its 
covering lattice (or line graph), the results here for bond 
percolation on the {7,3} lattice can be mapped to site 
percolation on its covering lattice, which is an interest- 
ing kind of hyperbolic kagome lattice. Covering lattices 
of the other lattices we considered here contain crossing 
bonds. 

Note added in proof: Just recently. Back [44 argued 
that the conjectured result Pu =0.5 for the EBT in Ref. 
[23] should be replaced by a lower bound of about 0.55, 
which is consistent with our results here. 
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Notes added after publication: This paper has 
been pubhshed in Physical Review E 85, 051141 (May 
29, 2012). We add the fohowing update and correction: 

1. Shortly after our paper was published, Delfosse and 
Zemor [45] posted a paper in which the rigorous upper 
bound for pi for the {5,5} lattice is reduced from 0.40 to 
0.38, closer to our approximate bound pf = 0.346. They 
also give upper bounds for self-dual lattices {m, m} for 
larger values of m. 

2. A correction in the printed version: In section III A, 



third paragraph, the second sentence should read: "We 
call this the duality point because on a truly dual lat- 
tice the crossing probability should also be 1/2. . . ." We 
have changed "occupation" to "crossing" in the present 
version. 

[45] Nicolas Delfosse and Gilles Zemor, "Upper Bounds 
on the Rate of Low Density Stabilizer Codes for the 
Quantum Erasure Channel," |arXiv:1205.7036| (May 31, 
2012). 



